Single polaron properties of the breathing-mode Hamiltonian 
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We investigate numerically various properties of the one-dimensional (ID) breathing-mode po- 
laron. We use an extension of a variational scheme to compute the energies and wave- functions of 
the two lowest-energy eigenstates for any momentum, as well as a scheme to compute directly the 
polaron's Green's function. We contrast these results with results for the ID Holstein polaron. In 
particular, we find that the crossover from a large to a small polaron is significantly sharper. Unlike 
for the Holstein model, at moderate and large couplings the breathing-mode polaron dispersion has 
non-monotonic dependence on the polaron momentum k. Neither of these aspects is revealed by a 
previous study based on the self-consistent Born approximation. 

PACS numbers: 71.38.-k, 72.10.Di, 63.20.Kr 



I. INTRODUCTION 

In a solid-state system, the interaction between a 
charge carrier and phonons (quantized lattice vibrations) 
leads to the formation of polarons. This mechanism is a 
key ingredient in the physics of the manganitesi and, 
possibly, of the cupratesi^'^ There are various model 
Hamiltonians describing the coupling of the particle and 
bosonic degrees of freedom. The asymptotic limits of 
weak or strong coupling can be investigated analytically 
using perturbation theory, however the intermediate- 
coupling regime generally requires numerical simulations. 
Recently, investigations of basic model Hamiltonians 
have progressed rapidly thanks to the development of 
efficient analytical and computational tools, and we are 
now able to begin to study more and more realistic mod- 
els. 

The simplest clcctron-phonon coupling is described by 
the Holstein Hamiltonian."* It is essentially the tight- 
binding model with an on-site energy proportional to the 



lattice displacement Xi = 



\/2MO 



H 



<ij> 



^r!^5^6,+5^n,X, (1) 



Here Ci is the annihilation operator for an electron at 
site i (since we only consider the single electron case, the 
spin is irrelevant and we drop its index in the following. 
We also set h — I), t is the hopping matrix, 



For the Einstein phonons, bi is the annihilation operator 
at site i, is the frequency, M is the atomic mass, and 
g is the electron-phonon coupling strength. The model 
has been widely studied numerically by Monte-Carlo 
calculations,—"— variational methods,'li^-22 and exact 
diagonalizationi^"— Analytic approximations have also 
progressed over the yearSi '^^'^^i'^^'^^ 

For some materials, a more appropriate model is pro- 
vided by the breathing-mode Hamiltonian. For example, 
consider a half-filled 2D copper-oxygen plane of a par- 
ent cuprate compound. Injection of an additional hole 
should fill an oxygen 2p orbital. Due to hybridization 
between oxygen 2p and copper 3d^2_y2 orbitals, the hole 



resides in fact in a so-called Zhang-Rice singlet (ZRS) 
with a binding energy proportional to — St^^, where tdp 
is the hopping between neighboring O and Cu orbitals. 
The dynamics of the ZRS can be described by an effec- 
tive one-band model with orbitals centered on the copper 
sub-lattice.**' If lattice vibrations are considered, the mo- 
tion of the lighter oxygen ions, which live on the bonds 
connecting Cu sites, are the most relevant. The hopping 
integral tdp and charge-transfer gap between Cu and O 
orbitals are now modulated as the oxygen moves closer or 
further from its neighboring Cu atom. Both the on-site 
energy and hopping integral are modulated in the effec- 
tive one-band model, but the former has been shown to 
be dominanl4ii^. The breathing-mode Hamiltonian de- 
scribes the physics of the linear modulation of on-site 
energy. 

While this breathing-mode model is motivated as a 2D 
model, in this work we investigate numerically only its ID 
version, relevant, e.g. for CuO chains. In ID we can in- 
vestigate accurately and efficiently not only ground-state 
(GS) properties, but also some excited state properties. 
For the Holstein model, it was found that polaron proper- 
ties are qualitatively similar in different dimensions,—"^ 
but with a sharper large-to-small polaron crossover in 
higher dimensions. We will show that a sharp crossover 
is already present in the ID model, and we expect less di- 
mensionality effects in the breathing-mode Hamiltonian. 

The ID breathing-mode Hamiltonian that we investi- 
gate here is described by: 
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V2Mn 



(2) 



The notation is the same as before, except that now the 
phonons live on an interlaced lattice. The difference be- 
tween the two models is more apparent in momentum 
space. The Holstein model has constant coupling to all 
phonon modes 



Vh - 
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whereas the breathing-mode model has a couphng 
strength that increases monotonically with increasing 
phonon momentum 



2ig 



^ sin |4_ Cfc {bl + 5_ 



Here N is the number of lattice sites, and becomes infi- 
nite in the thermodynamic limit. The momenta fc, q are 
restricted to the first Brillouin zone (— tt, tt] (we take the 
lattice constant a = 1). 

While numerical and analytical studies of the Hol- 
stein polaron abound, there is much less known about 
models with g{q) coupling. In particular, there is 
no detailed numerical study of the ID breathing-mode 
model, apart from an exact diagonalization of a sim- 
plified t-J-breathing-mode model in restricted basis^, 
and an investigation based on the self-consistent Born 
approximation^**, which is known to become inaccurate 
for intermediate and strong couplings. In this work, 
we study numerically various low-energy properties and 
the spectral function of the single polaron in the ID 
breathing-mode Hamiltonian. The results are compared 
with the relevant results for the single Holstein ID po- 
laron, allowing us to contrast the behavior of the po- 
larons in the two models. The article is organized as 
follows: in Sec. |TT]we review relevant asymptotic results, 
and describe the numerical methods we use to calculate 
low-energy properties and the spectral functions for both 
models. In Sec. IIIII we present our results, and in Sec. 
IIVI we present our conclusions. 



II. 



METHODOLOGY 



A. Strong-coupling perturbation results 

Perturbational results for the strong-coupling limit 
g ^ t provide a good intuitive picture of the problem 
even for the intermediate coupling regime. In the absence 
of hopping, t = 0, both Hamiltonians can be diagonalized 
by the Lang-Firsov transformation^ 



O = e^'Oe 



(3) 
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respectively, the diagonal forms of the Hamiltonians are, 
in terms of the original (undressed) operators: 
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(4) 
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\t{ni-nt+i) (5) 



i i 

where the kinetic energies are: 

3^ g(-''t+i + i.t) g(-ti+i+i,i) 



Q 

-te 2A7 
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For a d-dimensional lattice, Tb is modified by i) extra 
creation and annihilation operators of phonons in the di- 
rection transverse to hopping, and ii) change of the -3 
factor in the exponent to —{z + l),z — 2d. The third 
term in Eq. (jj]) and in Eq. ^ signifies that the mere 
presence of an electron would induce a lattice deforma- 
tion, leading to the formation of a polaron to lower the 
system's energy. For a single polaron, the lattice defor- 
mation energy is proportional to the number of nearest 
phonon sites (one for the Holstein model and z for the 
breathing- mode model). For t — 0, the ground-state en- 
ergy is degenerate over momentum-space: 



4°)(fc) = - 
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Each model has three energy scales, therefore the pa- 
rameter space can be characterized by two dimensionless 
ratios. It is natural to define the dimensionless (effective) 
coupling as the ratio of the lattice deformation energy to 
the free-electron ground-state energy —zt: 



A, 
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1 



2MV? zt 
1 



(6) 
(7) 



where z is also the number of nearest neighbors in the 
electron sublattice. It should be noted that since ~ 
l/V^M, the A's do not depend on the ion mass, M. A 
has been shown to be a good parameter to describe the 
large-to-small polaron crossover in the Holstein model. 
It will be shown in later sections that the definition also 
works well for the breathing-model model. The other 
parameter is the adiabatic ratio which appears naturally 
from the perturbation in t 



zt 



(8) 



Using standard perturbation theory,''^ the first-order 
corrections to the energy of the lowest state of momentum 
fc are: 



-2te-"^" cos(fc), 
,2te-i"^« cos(fc). 



(9) 
(10) 
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showing that the polaron bandwidth is exponentially sup- 
pressed in the strong coupling limit. As is well known, 
this is due to the many-phonons clouds created on the 
electron site (Holstein) or on the two phonon sites brack- 
eting the electron site (breathing- mode model). As the 
polaron moves from one site to the next, the overlaps be- 
tween the corresponding clouds become vanishingly small 
and therefore t^s — > 0. To first order in t, the suppres- 
sion is stronger for the breathing-mode model simply be- 
cause the overlap integral involves phonon clouds on 3 
sites instead of just 2, as in the case for Holstein. The 
second-order corrections are: 
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n 
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" fH{k,a,XH), 
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The functions / can be written in the form 



with 



fH,B = Ah,b + Bh,b cos(2fe) 



Ah = Ei{2aXH) - \n{2aXH) ~ 7 
Bh = Ei{aXH) — In(aAir) — 7 



and 



Ab — Ei{3aXB) ~ bi(3aXB) — J 



,aXf 



Bb = Ei{2aXB)-2Ei{aXB)+lni^)+l. 

Here, 7 is the Euler-Mascheroni constant and Ei{x) is 
the exponential integral with the series expansion 
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The result can be further simplified in the limit 

aXu.aXB > 1 using hm^_>oo I]^^! fiTT ~ ^- This 
leads to the simplified expressions: 
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^" cos(2fc) 



■ cos(2fc) 



(11) 
(12) 



Thus, the breathing-mode model's ground state en- 
ergy is slightly higher for any finite t. For the Holstein 
model, the dispersion is monotonic, since the second or- 
der cos(2fc) contribution is suppressed more strongly than 
the first order cos(fca) contribution. However, a quick 
comparison between Eq. (fTU]) and (|12p shows that in the 
breathing-mode model, the second order cos(2fc) contri- 
bution becomes dominant at large enough coupling. As a 
result, at strong couplings we expect the breathing- mode 
polaron energy to exhibit a maximum at a finite fc < tt, 
and then to fold back down. 



B. Matrix computation 

The computation method we use is a direct general- 
ization of the method introduced by Bonca et al. for the 
Holstein model, in Ref. ,171 . This approach requires sparse 
matrix computations to solve the problem. Although ex- 
pensive cpu and memory resources are required for this 
type of method, it gives us a systematic way to compute 
excited state properties, which would be more difficult to 
achieve using other numerical methods. 

The idea is to use a suitable basis in which to represent 
the Hamiltonian as a sparse matrix. For the Holstein 
model, this basis contains states of the general formii 



1 




(13) 



where K is the total momentum and S denotes a par- 
ticular phonon configuration, with sets of n„i phonons 
located at a distance m away from the electron. For 
the Holstein model m are integers, since the phonons are 
located on the same lattice as the electrons. The gen- 
eralization for the breathing-mode is simple: here m are 
half-integers, since here the phonons live on the interlaced 
sublattice. All states in either basis can be obtained by 
repeatedly applying the Hamiltonian to the free electron 
state which has all ri„i = 0. The possible matrix el- 
ements are —te^^^, fin, and ±-^=£=-^71, where n are 
integers related to the numbers of phonons. 

Since the Hilbert space of the problem is infinite, this 
basis must be truncated for computation. The original 
cut-off scheme in Ref. [l3 was optimized for computa- 
tion of ground-state properties of the Holstein model, by 
restricting the number of matrix elements between any 
included state and the free-electron state. Getting the 
higher-energy states is more involved, as it is evident 
from Eq. ([3]) that each state \(j>)LF in the Lang-Firsov 
basis correspond to a state 



)LF 



(14) 



in real space. With our choice for the S operators, this 
reverse transformation induces a phonon coherent state 
structure at the electron site (Holstein), respectively 
the two bracketing phonon sites (breathing-mode). The 
phonon statistics of the coherent state obeys the Pois- 
son distribution. In the anti-adiabatic regime {zt > fl), 
the splitting due to the hopping (off-diagonal hopping 
matrix elements) is significant compared to the diago- 
nal matrix elements proportional to Q. The underlying 
Lang-Firsov structure needs to be modeled by the hop- 
ping of the electron away from the coherent state struc- 
ture created by the e~'^ operator. To capture these char- 
acteristics, the basis is divided into subspaces with fixed 
numbers of phonons. Each subspace is enlarged by the 
addition of states with phonons further and further away 
from the electron site (increase of maximum value of m in 
Eq. (I13p ). until convergence is reached. This procedure 
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TABLE I: Most-relaxed cut-off condition for the ID 
breathing-mode polaron Green's function computation. 



Subspace's # phonons 


1 max 


# States 


1-11 


22.5 - (# phonons) 


17053356 


12-13 


10.5 


16871582 


14-15 


9.5 


28274774 


16-17 


8.5 


33423071 


18-20 


7.5 


41757650 


21-30 


5.5 


42628080 


31-40 


4.5 


38004428 


41-50 


3.5 


12857573 



allows for efficient generation of all basis states required 
to model the higher-energy states. 

Matrices of dimension up to 10^ — 10^ were needed to 
compute the two lowest-energy states accurately. These 
two states were calculated numerically using the Lanczos 
method with QR shift, ^^--^^ which works efficiently for the 
low energy bound states. 

The number of bytes required to store a n x n sparse 
complex matrix is roughly n((16 -I- 4)m + 4), where m is 
the number of matrix elements per row. The number of 
bytes required to store a n-vector is 16n. Therefore, an 
ordinary workstation can deal with n ^ 10^, sufficient for 
our low-energy states computation. Cluster paralleliza- 
tion provides decent speed-up up to n ~ 5 x lO'', above 
which communication costs proved to be too high due to 
the large matrix bandwidth, even after re-ordering. For 
the larger n values used in the Green's function calcula- 
tion (see below), SMP machines with high memory-to- 
cpu ratio were used for efficient computation. 

Table U lists the most- relaxed cut-off condition we used 
to calculate the Green's function (see next section) for 
the ID breathing-mode model. Because the queue time 
is roughly independent of memory requirements, but is 
longer than the computation time on the SMP machines, 
we relaxed the cut-off condition rather roughly until con- 
vergence was observed. As a result, the size of these 
matrices is certainly much larger than it has to be. 



C. Green's function computation 

Computation of higher-energy properties requires 
much larger matrices. The memory and fflops needed for 
such computations are formidable, especially to exten- 
sively investigate the multi-dimensional parameter space 
{X,a,K). Furthermore, one characteristic of the single 
polaron problem is a continuum of states starting at one 
phonon quantum above the K — groundstate. Lanczos- 
type methods typically are problematic in dealing with 
bands of eigenvalues with small separation. Therefore, in 
order to study higher energy states, we calculate directly 



the Green's functionis^ 

G(^,fc) = (fc| ^—\k), (15) 

where \k) = c^|0). This can be written as the solution of 
a linear system of equations: 

G{uj,k) = {k\y) 
{cu - H + ir])\y) = \k) 

One can iteratively tri-diagonalize H by the vanilla Lanc- 
zos process;^ 

H = QTQ^ 

iLo + iv-T)Q^\y) ^Q^\k) 

If the RHS is of the form [100...]^, Crammer's rule can be 
used to express G — {k\QQ^y) as a continuous fraction 
in terms of the matrix elements of the tri-diagonal matrix 
(uj + ir] — T). In particular, this condition is achieved by 
picking the initial Lanczos vector to be |fc). This method 
is efficient because it does not require the complete solu- 
tion of the linear system nor of the eigenvalue problem. 
It is well known that this type of iterative process suffers 
from numerical instability, which leads to the loss of or- 
thogonality in Q and incorrect eigenvalue multiplicity in 
y 50 -y^g perform the vanilla Lanczos tri-diagonalization 
and re-orthogonalize each states in Q against the starting 
vector |fc) to validate the continuous fraction expansion. 
Then, numerical errors may come from the fact that T 
may have the wrong eigenvalue multiplicity. However, 
this will not affect the location of poles in the spectral 
weight, i.e. the eigenenergies are accurate. 

III. RESULTS 
A. Low-energy states 

The ground state energy Eqs and quasiparticle (qp) 

weight Zq = |(<I>Gs|Cfc=olO)l^' where \^gs) is the ground- 
state eigenfunction, are shown in Fig. [1] for the ID 
breathing-mode and Holstein models. For a fixed value of 
a, we see the expected crossover from a large polaron (at 
weak coupling \b,h) to a small polaron (at strong cou- 
pling Xb,h), signaled by the collapse of the qp weight. 

The ground state energy of both models decreases 
monotonically with increasing coupling, but that of the 
Holstein polaron is lower. This is in agreement with 
the second order strong-coupling perturbation results in 
Eq. pT|) and p^ . Unlike the rather gradual decrease in 
the quasiparticle weight of the ID Holstein polaron, the 
ID breathing-mode polaron shows a large Zq at weak 
couplings, followed by a much sharper collapse at the 
crossover near As ~ 1. The reason for the enhanced 
Zq at weak couplings is straightforward to understand. 
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FIG. 1: (color online) (a) GS energy and (b) GS qp weight 
as a function of the corresponding dimensionless coupling pa- 
rameter. The dashed line corresponds to the Holstein model, 
while the breathing-mode results are shown by circles (line is 
a guide to the eye). Parameters are t = 2,Q — 1. 




FIG. 2: (color online) (a) Energy of the first excited K = 
state, measured from Eqs, and (b) its qp weight. The dashed 
line corresponds to the Holstein model, while the breathing- 
mode results are shown by circles (line is a guide to the eye) . 
Parameters are t = 2,Q — 1. 



Here, the wave-function is well described by a superpo- 
sition of the free electron and electron-plus-one-phonon 
states. Given the conservation of the total polaron mo- 
mentum K — = k + q and the large electron band- 
width t, states with high electron and phonon momenta 
have high energies and thus contribute little to the weak- 
coupling polaron ground-state. On the other hand, the 
coupling g(q) ^ sm(q/2) to the low-energy states with 
low electron and phonon momenta is very small for the 
breathing-mode model. This explains the slower transfer 
of spectral weight at weak coupling for breathing-mode 
versus the Holstein polaron. 

The energy (measured from Eos) and qp weight of 
the first excited K = state are shown in Fig. [21 For 
both models, at weak-coupling this state is precisely at 
O above the ground-state energy, at the lower edge of 
the polaron+onc-phonon continuum. As the coupling in- 
creases above a critical value, a second bound state gets 
pushed below the continuum. This second bound state is 
absent in SCBA calculations.^'* The separation between 
the two lowest energy states now first decreases and then 
increases back towards n as Xh,b oo. This behavior 
is well-known for the Holstein polaroniii The breathing- 
mode polaron shows the same qualitative behavior. Note 
that below the critical coupling, the computed energy of 
the first excited state is slightly larger than f2. The rea- 
son is a systematic error that can be reduced by increas- 
ing the number of one-phonon basis states, in order to 
better simulate the delocalized phonon that appears in 
this state. The qp weight of the first excited state is zero 
below the critical coupling, due to the crossing between 
on-site and off-site phonon statesJi 

The nature of the these states is revealed by checking 
the locality of the phonon cloud. We define the projection 



operator 

P{K)^ ^ \S,K){S,K\ 

Slocal 

where the summation is over all states with rriH = 
and niB — ±0.5 in Eq. (fTS]) . Comparison with Eq. p4|) 
shows that this operator selects only basis states with 
phonons only on the electron site (Holstein) and only 
on the two phonon sites bracketing the electron site 
(breathing- mode). Fig. [3] shows the expectation value of 
this operator for the two lowest eigenstates of both mod- 
els. For both ground-states (P(0)) ~ 1, indicating that 




nl , 1 , 1 , 1 , 1 

0.5 1 1.5 2 

\'\ 

FIG. 3: (color online) P(0) for t = 2, = 1 for the breathing- 
mode (red symbols) and Holstein (black line) models, respec- 
tively. The solid and dashed lines correspond to ground state 
and first excited state, respectively. 
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TABLE II: an_,n^ vs. n+,n_ for the ground state 



TABLE III: Q„ 



for the first-excited state 





t = 2 f7 


= 1 .9 = A 


= 




n+n_ 





1 


2 


3 





1.0000 


ooon 


nnoo 


onon 


1 


0.0000 


0.0000 


0.0000 


0.0000 


2 


0.0000 


0.0000 


0.0000 


0.0000 


3 


0.0000 


0.0000 


0.0000 


0.0000 






1 g = 1.5 A 


= 1.05 




n+n_ 





1 


2 


3 





0.9150 


0.0584 


0.1372 


-0.0477 


1 


-0.0584 


-0.1681 


0.0520 


-0.0536 


2 


0.1372 


-0.0520 


0.0549 


-0.3300 


3 


0.0477 


-0.0536 


0.0330 


-0.0255 




t = '2il = l 


g = 1.964 A 


= 1.964 




n+n- 





1 


2 


3 





0.9876 


-0.0509 


0.0082 


-0.0025 


1 


0.0509 


-0.0100 


0.0034 


-0.0020 


2 


0.0082 


-0.0034 


0.0022 


-0.0019 


3 


0.0025 


-0.0020 


0.0019 


-0.0019 



here most phonons are neariest to the electron. How- 
ever, at weak coupling the first excited state (which is 
here the band-edge of the polaron -I- free phonon contin- 
uum) has (-P(O)) 0, precisely because the free phonon 
can be anywhere in the system. When the second bound 
states form, (P(0)) becomes large, showing that phonons 
in these states are primarily localized near the electronjii 
While there appears to be a crossing between the ground- 
state and first-excited state values, we emphasize that 
P{K) measures the locality of the phonon cloud, not its 
structure. 

For the breathing-mode model, these results sug- 
gest the possibility to describe them using the on-site 
coherent-state structure. That is, a Lang-Firsov state 
with n_ and n+ number of phonons excited to the left 
and right of the electron site, mapped to real space by 
Eq. ([Ti)l . We note that we are no longer in the strong 
coupling regime and the transformation cannot be deter- 
mined by g and fl alone, therefore we seek an effective 
transformation with 

are projected into such 



The computed eigenstates 
structure an_,n+ by 




Table [TT] and IIIII show the results of such projec- 
tions for the ground state and for the first excited 
bound state. It is clear that they can be rather well 





t = 2 n = 


1 g = 1.5 A 


= 1.05 









1 


2 


3 





n 0228 


-0.6515 


0.0152 


0.1151 


1 


0.6515 


0.0062 


0.1700 


-0.0473 


2 


0.0152 


-0.1700 


0.0505 


-0.0456 


3 


0.1151 


-0.0473 


0.0456 


-0.0220 




t = 'lil^\ 


g = 1.964 A 


= 1.964 




n+n_ 





1 


2 


3 





-0.0859 


-0.6675 


0.0630 


-0.0164 


1 


0.6675 


-0.0864 


0.0246 


-0.0120 


2 


0.0630 


-0.0246 


0.0134 


-0.0099 


3 


0.0164 


-0.0120 


0.0099 


-0.0089 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 
k/K k/n 

FIG. 4: (color online) The fc-dependent (a) energy and (b) qp 
weight of the GS (circles) and first bound state (square) for 
the ID breathing mode model, for t = 2, n = 1, As = 1.07. 



e-^iW)^J2ie^'''cl\0), re- 



described as \GS) 

spectively |l)bound - e-^(*'f'"'^)^ e^^'cj (e'«5|_ ^ - 
6|_^i)|0) for some phase 9{t, g,n, K) needed to sat- 



-i9 



isfy time-reversal symmetry. These states no longer have 
definite parity symmetry like those of the Holstein model. 
The symmetry is broken by the anti-symmetric coupling 
term in the model. If the free-electron component is non- 
zero for an eigenstate, its components with odd/even 
number of phonons should have odd/even parity. For 
increasing momentum, this description is valid as long as 
the excited state remains bound, with energy less than 
Egs + ^■ 

Fig. |4] shows momentum dependent energy and qp 
weight for the two lowest eigenstates of the breathing- 
mode model, for an intermediate coupling strength Xb = 
1.07. The polaron band has a maximum at fc ~ 0.47r, in 
qualitative agreement with the strong-coupling pertur- 
bation theory results. This behavior is not captured by 
SCBA which is only accurate for low coupling strength4^ 
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FIG. 5: (color online) Ratio of effective polaron mass to that 
of the free electron. Circles and squares show breathing-mode 
results for GS and first bound state, respectively. The other 
lines correspond to the Holstein model GS (full) and second 
bound state (dashed). Parameters are i = 2, n = 1. 

For the Holstein polaron, the polaron dispersion is a 
monotonic function of momentum. Even though the 
qp weights remain moderately high at zero momentum, 
the weights collapse towards zero with increasing momen- 
tum, similar to the well-known Holstein case. This is due 
to the fact that at large total momentum, the significant 
contribution to the eigenstate comes from states with at 
least one or more phonons. The free electron state has a 
large energy for large momentum, and contributes very 
little to the lowest energy eigenstates, so indeed Z — > 0. 

The effective masses for the two lowest eigenstates of 
both models are shown in Fig. O as a function of Xh,b- 
These were calculated from the second derivative of the 
energy at momentum K — 0. For the GS of both mod- 
els, the effective mass increases monotonically with Xh,b- 
At weak couplings, the breathing-mode polaron is lighter 
than the Holstein polaron. As already discussed, this is 
due to the vanishingly weak coupling to low-momentum 
phonons. At strong coupling, however, the effective mass 
is larger for the breathing-mode polaron. This is in agree- 
ment with predictions of the strong-coupling perturba- 
tion theory, and results from the fact that the hopping 
of a breathing-mode polaron involves phonon clouds on 
2z — 1 phonon sites, whereas hopping of a Holstein po- 
laron involves phonons at only two sites. 

The effective mass of the first excited state can only be 
defined once this state has split-off from the continuum. 
It has non-monotonic behavior, first decreasing and then 
increasing with increasing Xh,b- This can be understood 
through the link of the effective mass and the qp weight. 
In terms of derivatives of the self-energy S(A;,a;), the ef- 
fective mass TO* is given by: 

TO* _ / as\ / TO a^EX"^ 

where derivatives are evaluated a,t K = and at the 
corresponding eigenenergy. The first term is linked to the 



qp weight, Z = (1 - f§)~^, so that m* ~ 1/Z. As shown 
in Fig. IHb), the qp weight of the first excited state has 
non-monotonic behavior, leading to the non-monotonic 
behavior of the effective mass. 

All the results shown so far were for a — 4. For higher 
a (lower fl and/or larger t), the difference between the 
two models can be grasped from Fig. \6\ Similar to the 
Holstein model, the large-to-small polaron transition oc- 
curs at lower A for increasing a^'^ At weak and mod- 
erate coupling, the qp weight and the effective mass (not 
shown) in the breathing-mode model are much less sen- 
sitive to an increase in a, than is the case for the Hol- 
stein polaron. This suggests that breathing-mode po- 
larons should be better charge carriers than the Holstein 
polarons, in this regime. 

B. Spectral Function 

The spectral function is proportional to the imaginary 
part of the Green's function: 

A{k,uj) = --lmG{k,uj) (16) 

TT 

In terms of single electron eigenstates and eigenfunction 
H\n) = En\n), we obtain the Lchmann representation: 

A(fc,c.) = 5]|(n|cl|0)p^(c.-£;„) 

n 

Of course, since we use a finite small 77 in numerical calcu- 
lations, the ^-functions are replaced Lorentzians of width 
T] [see Eq. (jisp j. We calculate the Green's functions as 
discussed in the previous section. 

Fig. [7] shows the spectral function for zero momentum 
as a function of energy. Results corresponding to four 
different coupling strengths Xb from near the crossover 
region are shown for the breathing-mode polaron. We 




FIG. 6: (color online) GS qp weights for the breathing-mode 
(circles) and Holstein (squares) models. Full symbols corre- 
spond to a = 8. For comparison purposes, the empty symbols 
show the a = 4 results of Fig. [Jb). = 1) 
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FIG. 7: The spectral function of the ID breathing-mode 
model, on linear (left) and logarithmic (right) scales. The 
curves are shifted for better viewing, and correspond (top to 
bottom) to As = 1.5; 1.25; 1.0 and 0.75. Vertical lines indi- 
cates Egs and Eos + K = 0,t = 2,n ^ l,ri ^ O.OOAfl. 



note that there is always a continuum starting at one 
phonon quantum above the ground state energy. This is 
more clearly visible in the right panel, where the spectral 
weight is shown on a logarithmic scale, and vertical lines 
mark the position of the ground-state energy -Egs, re- 
spectively of Eos + ^- As the couphng A_b increases, we 
see the appearance of the second bound state below the 
continuum. We only find at most one extra bound state 
in this energy range for all coupling strengths. As Ab 
increases, the spectral weight of the first continuum de- 
creases dramatically. Other bound states form above it, 
followed by higher energy continua whose weight is also 
systematically suppressed. This is qualitatively similar 
to the behavior exhibited by Holstein polaron.'^^ 

Fig. [5] illustrates the momentum dependence of the 
breathing-mode polaron's spectral weight. The results 
correspond to a coupling above the critical value, where 
there is a second bound state. The majority of the spec- 
tral weight is transferred to much higher energies as the 
momentum increases, and a broad feature develops at 
roughly the position of the free-electron energy for that 
momentum. This spectral weight transfer is also qualita- 
tively similar to what is observed for Holstein polarons. 
Our results have a high-enough resolution to clearly show 
the continuum at Egs + ^ for all values of K. This is 
part of the kink-like structure reported in Ref. |3 The 
logarithmic plot clearly reveals a non-monotonic disper- 
sion of the ground-state like in Fig. [H characteristic for 
the breathing-mode polaron. 

Fig. [5] shows only one peak located between the 
ground-state and the polaron+one-phonon continuum, 
even though in fact we believe that there are more than 
one eigenstates within this region. We found, from eigen- 
value computation, additional energy states below the 
continuum; however, computation of exact energy val- 
ues requires prohibitively long compution time due to 
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FIG. 8: A{k,uj) vs uj, for K/tt^O; 0.25; 0.5; 0.75 and 1 (top 
to bottom) for intermediate coupling t = 2,i} = 1,Xb ~ 1.25. 
rj/4: = AE = O.OOm. The height of the spectral weight is 
plotted in linear scale on the left and logarithmic scale on the 
right. The two vertical lines indicate Egs and Ecs + ^■ 



the clustering of eigenvalues. By observing the conver- 
gence behavior due to increasing basis size, we can con- 
clude that additional bound states do exist below the 
continuum. The lack of their contribution to the spec- 
tral function can be understood by the fact that the sin- 
gle particle Green's function only contains information 
about eigenstates with finite qp weight, KfAlc^lO)!^ > 0, 
see Eq. ([15]). These eigenstates must have components 
corresponding to some Lang-Firsov eigenstate with no 
off-site phonons (Eq. (fT4|) ). Also, the wave- function of 
these states must have a peculiar space inversion sym- 
metry: S-symmetric for all even-phonon-number compo- 
nents and P-symmetric for all odd-phonon-number com- 
ponents. The ground-state always satisfies this require- 
ment, but above a critical coupling, only one other state 
below the phonon threshold satisfies this requirement. 



IV. CONCLUSIONS 

In summary, we have reported here the first accurate 
numerical study of the ID breathing-mode polaron. A 
previous study^ based on the self-consistent Born ap- 
proximation proves to be inadequate to describe correctly 
the behavior for medium and strong couplings, as ex- 
pected on general grounds. 

Comparison with the Holstein model results, which 
correspond to a coupling g(q) — const., reveals some of 
the similarities and differences of the two models. The 
breathing-mode polaron is much more robust (has a much 
larger qp weight, and less variation with parameters) at 
weak couplings. This is a direct consequence of the fact 
that coupling to low-momentum phonons, which is rele- 
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vant here, becomes vanishingly small g{q) ^ sin(q/2) 
0. Similar behavior is expected for any other g{q) model if 
limq^o 9{q) 0. On the other hand, at strong couplings 
the breathing-mode polaron is much heavier and has a 
lower qp weight than the Holstcin polaron. This also re- 
sults from strong-coupling perturbation results, and is 
due to the fact that in order to move from site i to site 
i -|- 1, a small breathing-mode polaron must (i) create a 
new polaron cloud at site i + 1 ; (ii) rearrange the polaron 
cloud at site * + 5 , so that its displacement is now point- 
ing towards site i + 1, not towards site i; and (iii) relax 
(remove) the phonon cloud at site i — ^- This results 
in a suppressed polaron kinetic energy, and an enhanced 
effective mass. Because of the larger Z at weak coupling, 
and the lower Z at strong couplings, the crossover from 
large to small polaron is much sharper for the breathing- 
mode polaron. Another interesting observation is that 
the polaron dispersion becomes non-monotonic with mo- 
mentum k for medium and large couplings. This can 
be understood in terms of strong-coupling perturbation 
theory, which shows that the second-order cos(2fc) cor- 
rection is larger than the first-order cos(fc) correction for 



large-enough couplings. 

Similarities with the Holstcin behavior regard the 
appearance of the polaron-f free phonon continuum at 
Eqs + and the appearance of a second bound state 
with finite qp weight for large-enough couplings. The con- 
vergence of numerics points to the existence of additional 
bound states, whose absence from the spectral function 
can be explained by symmetry or missing free-electron 
componenets in the wavefunction; however, this issue is 
not fully settled. Also, the importance of such states for 
the physical properties is not known. The general aspect 
of the higher energy spectral weight at strong couplings, 
as a succession of bound states with large spectral weight 
and continua with less and less spectral weight, is also 
reminiscent of the Holstein polaron results. 
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